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This  report  falls  naturally  into  two  parts  and  is  written  in  that  form. 
The  first  deals  with  fracture,  the  second  with  matters  relating  to  metallic 
grains. 

PART  I  FRACTURE 

The  object  here  has  been  to  develop  mathematical  models,  involving 
singular  integral  equations,  for  cracks  in  the  presence  of  plastic  deforma¬ 
tion  and  to  use  these  models  so  as  to  predict  the  stress  necessary  for 
crack  propagation. 

It  is  well  known  that  the  processes  of  plastic  deformation  in  the 
neighborhood  of  a  crack  impedes  its  advance.  It  has  not  been  the  role  of 
fracture  mechanics  to  provide  an  understanding  of  why  this  should  be,  but 
rather  to  give  a  workable  system  of  prediction  of  fracture  behavior  from 
certain  particular  data.  To  date,  considerable  success  has  been  achieved 
by  following  this  route.  However,  this  success  has  not  been  uniform  and 
present  achievement  falls  short  of  the  desired  goals. 

This  situation  has  prompted  more  fundamental  work  of  which  that 
supported  here  represents  a  further  stage  of  advance. 

An  essential  deficiency  of  fracture  mechanics  is  that  the  plastic  zone  at 
the  tip  of  a  crack  in  a  plastically  deformable  material  is  treated  simply 
as  a  region  in  which  elasticity  is  non-linear.  This  gives  rise  to  a  major 
conceptual  difficulty,  the  so-called  "Rice  paradox,"  namely:  that  as  a 
consequence  of  the  need  to  consider  an  elastic-plastic  continuum  there  is 
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no  force  available  to  provide  the  surface  energy  of  a  crack  as  it  advances. 

Thus,  fracture  mechanics  gives  an  imperfect  representation  of  the 
physics  of  fracture.  Again,  we  note  that  it  does  not  take  accoint  of  the 
role  of  dislocations  in  achieving  plastic  deformation,  nor  does  it  allow 
for  the  interactions  which  occur  between  these  dislocations  nor  for  that 
between  them  and  the  crack. 

An  immediate  consequence  of  the  representation  of  a  plastic  zone  in 
terms  of  dislocations  is  the  realization  that  there  must  be  a  finite  region 
surrounding  the  crack  tip  which  is  free  of  dislocations.  This  realization 
allows  the  identification  of  the  source  of  the  Rice  paradox  and  it  is  seen 
that  it  is  vital  towards  a  thoroughgoing  treatment  of  fracture  to  include 
such  a  dislocation-free  zone  in  the  analysis. 

As  a  first  exercise  in  this  direction,  Chang  and  Ohr  analysed  the 
behavior  of  a  finite  planer  two-dimensional  crack  with  two  plastic  zones  in 
the  plane  of  the  crack.  Thus,  the  crack  was  supposed  to  lie  in  the  plane 
y  =  0  in  the  region  -c<x  <c  with  plastic  zones  in  the  regions  b<|x|<a 
with  a>b>c. 

Representing  equilibrium  in  this  situation  in  terms  cf  a  singular  inge- 
gral  equation  ^  Chang  and  Ohr  were  able,  in  particular,  to  find  the 
dislocation  distribution  function  and  thence  the  force  for  crack  extension. 

The  geometry  adopted  by  Chang  and  Ohr  should  be  appropriate  for  frac¬ 
ture  in  modes  II  or  III,  but  these  are  not  generally  situations  of  technical 
interest.  This  deficiency  has  in  part  been  remedied  in  the  course  of  the 
present  work. 

The  essential  difficulty  encountered  in  extending  this  approach  to  more 
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realistic  geometries,  namely  those  in  which  crack  dislocation  are  skew 
one  to  the  other,  is  that  the  governing  integral  equations  become  much 
more  complex.  This  difficulty  has  necessitated  the  development  of  a 
method  of  solution.  This  method  is,  of  course,  new  and  its  development 
is  a  major  component  of  the  results  of  this  research.  It  is  detailed 
below. 

Skewed  D.F.Z.  Cracks  in  Modes  One  and  Three 

The  complications  in  the  analysis  of  D.F.Z.  cracks  are  such  that  when 

121 

solutions  in  closed  form  are  available1  1  the  results  can  be  expressed  only 

in  terms  of  higher  transcendental  functions  whose  significance  becomes 

transparent  only  after  approximations  are  made.  Accordingly,  it  is 

germane  to  consider  the  use  of  treatments,  which  while  approximate,  are 

nevertheless  quite  accurate  and  allow  the  analysis  to  be  carried  through 

and  the  final  result  expressed  in  terms  of  elementary  functions.  This 

[41 

approach  has  been  used  previously  for  the  simplest  case1  J  where  the 
crack  and  its  associated  plastic  zones  are  col  inear.  It  is  here  used  to 
treat  the  more  complicated  and  physically  more  realistic  situations  in 
which  the  plastic  zones  are  skew  to  the  crack  plane. 

The  particular  method  employed  is  also  applicable  to  the  colinear  case 
where  it  is  most  simple.  It  is  as  set  out  below: 

We  consider  a  planar  two-dimensional  crack  -c<x<c  in  a  uniform  iso¬ 
tropic  medium  uniformly  loaded  so  as  to  provide  a  shear  stress  a  in  the 
plane  of  the  crack.  Plastic  zones  are  supposed  to  lie  in  the  region  b<x<a 
and  to  be  represented  by  continuous  distributions  of  dislocations  formed  in 
the  presence  of  an  applied  stress  in  the  four  straight  lines  which  make 
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angles  +y  and  +(7r-y)  to  the  plane  of  the  crack.  It  is  supposed  that  the 
operative  applied  stress  in  these  regions  is  o-r,  (r>o). 

We  calculate  and  solve, 

(a)  the  dislocation  distribution  g(x),  which  vanishes  at  |x|=b,  |x|=a  and 
which  is  developed  in  the  regions  b<|x|<a  due  to  a  dislocation  dipole 
with  elements  of  strength  +f  (or)  located  at  x  =  +,  cr<c; 

(b)  the  magnitude  of  a  uniform  stress  Ar(cr)  which  must  be  supposed 
present  in  order  that  the  dislocation  distribution  exist; 

(c)  the  stress  developed  in  the  region  -c<x<c  by  the  distribution  g(x); 

(d)  the  resulting  singular  integral  equation  for  the  dislocation  distri¬ 
bution  representing  the  crack  with  the  constraint  that  the  implied  sum 
over  the  elements  f(a)  gives  ZAr(ar)  =  r. 

In  cases  where  the  distributions  representing  the  plastic  zones  are 
skew  to  the  crack  length  the  interactions  between  dislocation  elements  are 
not  so  simply  expressed  and  this  leads  to  some  complications  and  inaccur¬ 
acy  in  the  analysis.  We  commence  with  an  examination  of  these  inter¬ 
actions.  Referring  to  Fig.  1 


e 


FIG.  1 

Crack  Tip  Geometry 
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we  are  interested  to  find  the  forces  acting  in  the  direction  OA  and  OC  on 
dislocations  magnitude  Xt,  X2  located  at  the  points  C  and  D  distant 
respectively  q  and  £  from  the  point  O. 

We  *.wO  cases,  the  dislocations  are  (A)  screw  and  (B)  edge. 

(A)  Screw  Dislocation 

In  this  case,  the  O-B  component  of  the  force  acting  on  the  dislocation 
at  D  is 

F  =  Ag  Xt  X2  cos  <p/X2q 

where  A  =  /A/2tt.  From  geometry  it  is  readily  shown  that 
© 

£  +  q  cos  y  (1) 

rf  +  ^  +2q(  cos  y 


F  = 


AgXiX2 


AgXj  X2 

r)  cos  y 

X' 

+ 

XtX2 

f  l 

_  (  1  -  cos  y  ) 

X* 

l  n +  f 

n  ~lq  +  ZV 

(2) 


>  AsX,X2  r 

I*  V 


COS 


+ 


LAX 


<n 


cos 

~TF 


jd 


► 

4 


with  a  fractional  error  which  is  greatest  when  q  =  £  and  2r)£(l-cos  y)/ 
(rj  +£)2  =1  which  then  underestimates  F  by  a  factor  which  is  25% 
when  y  =  60°. 
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(6)  Edge  Dislocations 

Supposing  the  Burgers  vector  to  be  parallel  to  the  lines  on  which  they 
are  placed,  the  force  which  acts  along  the  line  of  grain  between  the 
dislocations  is  given  by 


F  = 


yuXiX2 


r\  2tt(  1  -i-' )  r\  ~  A2 


e  XjA2  cos 


and  the  contribution  to  the  force  for  glide  along  the  line  OB  is 


c-  ^  _  ^e^1^2  cos  y  cos 
Fq  cos  <p  =  — p —  ^ 


Additionally,  there  is  a  component  of  force  from  Fg 


r-  ^  _  ^e^1^2  sin  A  sin  <t> 
Ffl  sin  <p  -  -YZT-- 


There  results  a  total  glide  force 

^e^1^2  cos  (y -<p) 

X2r} 


On  the  same  basis  as  that  used  for  screw  dislocations  we  find  a  glide 
force,  e.g.  (1) 

-  {“»  '  n  +  sln2y  ^  0!  } 

=  ^  fcV+Vln;jf)  -  TFT?  0’  < r  ■ -  «*»  -  s,n’r,  } 


or 
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AgXjX2 
-  ' 


{  {  cos  jn^-cosjd  +  sin2Y 


} 


AgA^j  |"  j  +  £  cos  y  -  C  cos2y  -  r>  sin2y  1 

7?  IT+T  Tn“+#  I 


or  better 
A_AjA2 


MeAl  2  r  cos2y  +  cos  y  £  (1  ‘  cos  y)  +  sin2y  f  £  11 

-  ~^r~  t n+T  Oi  +  tf2  tff+Url/ 


_  AgA^ 

“X2 


f  cos2y  +  £  (cos  y  -  cos2y)  1 

ln  +  e  (0  +  ^)2  / 


We  are  now  in  a  position  to  proceed  with  the  sequence  of  calculations 
SDecified  above  which  leads  to  a  determination  of  the  dislocation  distri- 

A 

bution  for  the  crack.  We  consider  the  case  of  mode  III  loading  and  so 
screw  dislocations. 


Distribution  Function  For  Mode  III  Loading 

Referred  to  ordinary  coordinates  the  force  acting  on  a  dislocation, 
Burgers  vector,  X2  at  a  point  x  =  t  due  to  a  dislocation  at  q  =  a  is  from 
(1): 


F  = 


X,  x2 

T7- 


{t 


a 


g(t-c) 
(t-a  )  2 


} 


Hence  to  facilitate  the  calculation  we  make  another  approximation.  We  see 
that 

0  (t-c)  -  0  t(t+c)  (t-c)  =  0  t  (t2-c2) 

“^d7-  It^I2  IcHFaP 

with  error  for  a  plastic  zone  size  A  -  small  compared  with  the  crack 


length  (i.e.,  A/c<l)  gives  an  overestimate  of  the  term  in  0  of  at  most  5% 
and  partially  offsets  the  underestimate  involved  in  the  previous  approxima¬ 
tion. 

On  this  basis  we  calculate  the  dislocation  distribution  due  to  a  pair  of 
dislocations  +  A  located  at  points  x  =  ±or  together  with  a  uniform  applied 
stress  <5r.  The  equation  of  equilibrium  for  this  distribution  is 


i  >  dx 

St  +  A_  A 

1 

-  i 

gM  x.t 

Jo 

SX*  l 

t  -  a 

t  +  a J 

The  distribution  is  chosen  to  vanish  at  each  of  the  end  points  x  =  +b,  +a. 

[1] 


Following  Head  and  Louat  L  J  we  have 

gM=  "A  |V-x2)  (*2-b2)  1  1  1 


/  JL-  + 
1  t-or 


D  [(a2-t2)  (t2-b2)]^ 

/ 3 1  (t2  -  c2)  +  6x  -  p  -  §  t  (t2  -  c2)  \  dt 
2c2  (t  -  or)2  t  +  a  2c7  (t  +  a)2  J  t-x 


which  rearranged,  becomes: 


eM  = '  A 


fa2-r)2)  (rj2-b2) 


i 


[(a2-t2)(t2-b2)]*  ^ 


£  3o 

tz-o 

+  St 

dt 

- 

t-x 

where  $  =  (1  -  cos  y)  /c2 


Here  integration  is  to  be  carried  out  over  the  union  of  segments  b<|t|<a 
and  is  easily  effected  by  the  use  of  the  appropriate  contour  and  the  calculus 
of  residues.  We  find: 
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(  ]  _  A  f(a2-x2)  (x2-b2)~|  ’  {  ap  +  p  — g  '  C  ^  ~  l  '  nj~. 

gM  -  L  J  l  °a  J  [(a2-a2)  (b2-a2)]*  Ca 
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provided  1 


x) 


1  _  f 6x  +  A  A  f  a  +  /S’  t2(t2  -  c2)  d  1  1  \  ts  dt  =  0; 

[(a2-t2)(t2-b2)]*  ^  ^ 

s  =  -  0,  1,2,3. 


This  is  satisfied  from  consideration  of  symmetry  for  s  =  0  and  s  =  2 
so  that  we  are  left  with  two  conditions  associated  with  cases  s  =  1  and  3. 
These  conditions  may  be  written  as: 


~T 


6v  a2  +  b2  +  crpk  +  a/5k(a2  +b2  +  4or2  -  2c2) 

/  kp  +  k/S  (5as  -  4c2  -a2  (c2  -  a2)  \  ,  s  =  3; 

3  1  (a'  -  a2]  J 

_  r  .  .  .  .  ^  1 

\  »  s  =  1, 


-  a3 


[(a2-a2)  (b2~a2)] 
a2  Sr  =  - 


[(a2-a2)  (b2-ar2)  j 


73 


f  kp  +  k/S 

r  \  N 

3a2  -  2c2  -  a2  (c2  -a2) 

l 

^  a *  -  a2  '  ' 

where 


k  =  2AA. 


~X~ 

We  see,  for  the  important  range  where  a/c  =  0  (1),  that  the  right  hand 
sides  of  these  equations  are  essentially  equal.  Accordingly,  we  find  that 


6z  *  -2kor  [p  +  P  (a2  +  b2  +  4a2  -  2c2)  ] 
(V  +  b2  “2a25 - 


(3) 


provided  a/c  ~  1,  We  shall  see  later  that  this  approximation  is  adequate 
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for  our  purposes. 

We  now  determine  the  stress,  in  the  plane  of  the  crack,  developed  by 
the  dislocations  of  the  distribution  g(x).  We  have,  using  an  approximation 
similar  to  that  used  before,  a  force  on  a  dislocation  X,  at  x,  directed  along 
the  crack  plane,  developed  from  another  dislocation  X2  located  at  some 
distance  t  up  the  skewed  slip  plane: 


F=  A_  X,X?  f  1  -  jS(t  -  c)  ]  ~  A_X,X2  f  1  -$t  (t2-c2)  )  (4) 

-3^-?!*—  {—t  'p  (x .,7*-  ) 

Then  using  (2)  and  (3)  we  have  at  a  point  x  along  the  crack  a  stress 

cr(x)  =  -  AAg  f  ap  +  d  a2(cr2  -  c2)  j  1 _ 

~X~  L  '  ((a2>a2)  (b2-cr2)]^ 


JD  0  '  1 


5  t(t2  -  c2)  1  dt 


where  integration  is  again  over  the  union  of  the  intervals  -a<t<-b,  b<t<a. 
Carrying  out  this  integration  we  have: 


<7(x,cr)  =  ttAA  f  ap  +  £  d  cr2(a2  -  c2)  \  1 _ 

35  *  [(a2-a2)  (b2-cr2)]^ 

I  1  +  [(a2-x2)  (b2-x2)]^  -  f(a2-g2)(b2-g2)]  ^  I 

|_  v*  -  ni  -I 


-  jB  d  /  x3  +  x  f  f(a2-x2)(b2-x2)]^  (c2-x2)  -  f(a2-or2)  (b^or2)]^  (c^cr2) 

icl  x2  -  a* 


Now  this  is  the  stress  resulting  from  a  single  dislocation  A  located  at 
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4 


i 


I 


i 


1 


i 


a.  The  total  stress  developed  at  a  point  u  on  the  real  axis  by  a  disloca¬ 
tion  distribution  f (or)  over  the  length  of  the  crack  is 


f  f  (a)  f  1  +  a(u,a)  ] 

t  \u=a  -rir) 


da 


Now  since  the  crack  constitutes  a  free  surface  this  stress  together  with 
the  applied  shear  is  to  vanish. 

Then  on  the  basis  that  the  distribution  is  unbounded  at  t  =  +c  we  have 

[4]: 


f(u)  = 


4  r_  fC  t  dx_  _  A  fc 

- il  [c2-x2]2  x-u  £  [c2-n2] 

ir2A[c2-u2]^  J  J 


i 


n-x 


f(a)  a(q,a)  da 


} 


Performing  the  first  integral  and  interchanging  the  order  in  the  second, 
this  becomes: 


f(u)  =  1 


7r[c2-U2] 


i 


cru 

“A 


1 

2 


rc 


{ 


f(a)  f  ap  +  d  a2(a2  -  c2) 

Ha 


} 


1 


[(a2-a2)  (b2-a2)J 


rC 

1  +  [(a2-o2)  (b2-n2)]  ^  -  [(a2-a2)  (b2-a2)]  ^ 

J-C  - 

•  f  d  f  +  x  !  ffa;-r)2)  (c*-n2)  -  (c*-a*) 

aa  \  1  q*  -  a 7 


4 


fcV)1  £u]dr> 
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We  now  introduce  the  approximation 

[{a2-/}2) (b2-rj2)]^  ~  ab-k^2  ,  k  =  ab  -  [(a2-b2) (b2-c2)]^ 


15) 


thereby  avoiding  the  complications  associated  with  elliptic  integrals.  The 
error  introduced  in  so  doing  is  estimated  (graphically)  to  be  less  than  5% 
and  to  be  such  that  a  reduction  of  k  by  such  an  amount  would  tend  to 
eliminate  that  error. 

Integrating  in  (5)  we  have 


f(u)  =  u 


rr[c2-y2] 


1 


rc 

1  !  f  (or) 

A 


|  orp  +  j9_d_  o2(ar2-c2)  j 


(6) 


k+1 
[(a2 


-f-i  ftc2  f  (3  -  k)  u2  +  3ab  ■  k  |  da 
-a2)  (b2-cr2)]7  <  FT7"; 


Then  remembering  that  we  will  be  concerned  only  with  the  behavior  of  f(u) 
as  u-k:  we  simplify  matters  by  writing 

(3-k)  u2  =  3-k 

F 

in  (6). 

On  this  basis  we  then  recognize  that  (5)  can  be  written  as 

f(u)  =  u  /  o  -  F(a,b,c)  \ 

t  1  ?’  i  \  IT  A  J 

Icz-UzJ  z  ^  J 

where  F(a,b»c)  can  be  evaluated  by  substituting  this  expression  for  f(u)  in 
(6).  We  get 


7T[c2-U2]: 


[a  -  i 
IS  Z 


1  q  f  a  -  F(a,b,c)  1 
5  rXI?  I  tFA  / 


_c[c2-c*l 


f  ap  +  d  q2(a2-c2)  1  k+1  -  ffc2  f  3-k  +  3ab  -  k  \  da 

\  J  [(a2-a2)  (b2-a2)]^  ^  ~sr  ) 


F(a,b,c)  i 


tt[c2-u2)]  ’ 


Thence  we  have: 


F(a,b,c)  =  <7  G(a,b,c,) 

7rAs  G(a,o,c}-Z 


where 

G(a,b,c) 

rC 

~  q  |  ap  +$  d  o2(a2-c2)  k+l-ftc2  /  3  +  ®b  "|  -2k  j 

_c[c2.a2ji  l  0a  J  [{a2-a2)  {b2_a2]*  L  l  ^  J  J 


Once  again,  to  facilitate  integration  we  use  the  same  approximation 
[(a2-x2)(b2-x2)]^  =  ab  -  kx2 


and  find 


G(a,b,c)  =  n  D  [ab] * 
k[ab-kc2]‘ 

where  D  = 

(k+1  -  £c2)  (3-k  +  3ab  -  k) 
c2 


Then  substituting  for  G(a,b,c)  we  hcve 


aP  [ab]^  {  p  +  0ab  } 

F(a,b,c)  =  jS  _ 

^  D[ab]  *  {p  +  /5ab}  -  k[ab-kc2]*  j 

Here  wa  have  neglected  terms  not  involving  the  factor  (ab-bc2)"^  which 
gives  a  large  contribution  when  the  dislocation  free  zone  is  relatively 
small,  i.e.,  when  [c(b-c)^  >>1. 

It  remains  to  relate  the  uniform  stress  r  which  acts  to  oppose  dis¬ 
location  motion  in  the  intervals  b>|xl<a  of  the  plastic  zones.  Using  (3), 
we  replace  A  by  f(x)6a  and  integrate  over  the  length  of  the  crack  to  obtain 


r  = 


-  4A 

IT 


[Cf(or)a 


J-c 


f  p  +  £(a2  +  b2  +  4a2  -  2c2) 
\  a2  +  V - ia2 - 


} 


da 


=  -  4AH 

1FT 


[c2-a2)  ^ 


da 

(a2+b2-2a2) 


where 

H  =  {  p  +  £(a2  +  b2  +  4a2  -  2c2)  ]  {  £  ‘  F-' (a’b’c)  ] 


Integrating  we  have 


x-n H 

"IT 


a2  ^  b2 
b2  +  a2  -  2c2 


(9) 


This  completes  the  solution.  This  is  not  exact  since  it  includes  a 
number  of  approximations.  However,  none  of  these  leads  individually  to 
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serious  error.  Furthermore,  they  appear  to  broadly  balance  in  over¬ 
estimating  in  some  cases  and  underestimating  in  others  so  that  the  net 
error  should  be  small,  i.,e.,  a  few  per  cent. 

The  quantity  important  in  considerations  of  fracture,  namely  the  crack 
extension  force  G  (the  square  of  k),  follows  immediately  from  (7)  through 
the  use  of  the  Bilby-Eshelby  formula 

G  =  ttA  A  lim  (u-c)  f(u) 

~T~  u -*=  • 

Using  this  result  together  with  (8)  and  (9),  G  can  be  enumerated  for 
cases  of  interest. 

The  Interaction  Of  Cracks  With  Cylindrical  Microvoids 

A  characteristic  of  fracture  when  significant  amounts  of  plastic 
deformation  occur  near  the  tips  of  a  crack  is  the  appearance  of  micro¬ 
voids.  These  microvoids  are  roughly  circular  in  cross  section  and 
centered  in  the  plane  of  the  crack.  Once  these  voids  are  formed  further 
crack  extension  appears  to  involve  the  coalescence  of  crack  and  void. 

Some  part  of  this  process  is  achieved  by  crack  motion.  It  is,  therefore, 
of  interest  to  examine  how  the  crack  extension  force  (G)  is  affected  by  the 
presence  of  the  microvoid. 

Towards  this  end,  we  consider  a  somewhat  simplified  two-dimensional 
representation  in  which  the  crack  is  confined  to  a  plane  and  the  microvoid 
is  a  cylinder.  Specifically,  we  suppose  the  crack  to  lie  in  the  plane  y  =  0, 
in  the  range  -c<x<c  and  to  be  unbounded  in  the  z-direction.  We  suppose  the 
crack  behaves  symmetrically  so  that  microvoids  are  located  at  each  end. 
These  are  taken  to  be  cylinders  of  radius  r  having  their  axes  parallel  to 
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the  z-axis  and  passing  through  the  points  a*  o,  and  -a,  o. 

In  its  full  generality  the  determination  of  G  for  this  geometry  is  a  very 
complicated  process  for  which  it  is  improbable  that  an  analytical  solution 
is  available.  However,  such  generality  is  not  essential  here  where  r<<c. 
Thus,  we  suppose  that  the  microvoid  behaves  in  the  stress  field  of  crack  as 
though  it  were  subjected  to  a  uniform  stress  aQ  having  the  value  of  a 
developed  at  the  position  of  the  center  of  the  void  when  the  crack  alone  is 
present.  The  approximation  so  involved  has  been  investigated  and  shown  to 
be  viable  [6]  for  the  case  where  the  void  is  replaced  by  a  thin  microcrack. 

This  behavior  can  be  traced  to  the  fundamental  feature  that  the  stress 
variation  over  the  length  of  a  microcrack  is  for  the  main  part  represen¬ 
table  by  odd  functions  of  distances  referred  to  the  center  of  the  micro¬ 
crack.  These  odd  functions  are  found  to  be  relatively  unimportant  in  deter¬ 
mining  the  far-field  stresses  developed  by  the  microcrack.  The  reason 
for  this  characteristic  can  be  traced  to  the  requirement  that  the  total 
dislocation  content  of  the  microcrack  be  zero.  A  result  of  this  constraint 
is  that  the  far-field  stresses  exhibit  quadrapolar  for  odd,  as  opposed  to 
dipolar  for  even,  functional  behavior.  Since  the  same  constraint  must 
apply  for  microvoids  we  can  expect  them  to  exhibit  similar  behavior. 

Again,  it  is  possible  to  make  a  simple  examination  of  the  effects  of  the 
stress  components  representable  by  even  functions.  Thus,  we  write  the 
stress  in  the  neighborhood  of  a  crack  under  a  uniform  stress  a  as: 

2  =  cm  =  cr(a+t) 

[r)2-c2]^  [(a+t)2-c2]^ 


(11) 


Thus,  we  suppose  the  void  to  be  centered  at  rj  =  a  and  that  the  variable  t 
measures  distances  from  that  point.  Expanding  this  expression  in  a  Taylor 
series  we  find  a  mean  stress  over  the  interval  -A<t<A: 


2  =  era  f  1  +  c2A2  +  3c2  (4a2  +  3c2)  A4  +  —  1 

l~^T  1  ^T“WF  J 

For  A<>  •  7(a-c)  and  c~a  this  series  is  rapidly  convergent  and  Z  nearly 
constant  and  equal  to 


I  =  era 

[a2-c2)]* 

which  is  the  stress  calculated  from  (1)  at  the  position  of  the  center  of  the 
void.  For  the  reasonable  case  where  A  =  r  =  a-c/2  the  actual  mean  value 
deviates  from  Z  by  only  about  3%.  Accordingly,  we  proceed  as  though 
the  microvoid  were  subjected  to  the  uniform  stress 

Now,  if  the  dislocation  distribution  which  represents  the  rate  of  change 
of  surface  displacements  over  the  crack  surface  is  f(x),  the  stress  normal 
to  the  plane  of  the  crack  at  a  point  t,o  due  to  the  distribution,  is  given  by 


-A 


(10) 


Here  A  =  fiA/27r(l-y),  p  is  the  shear  modulus  for  the  material  assumed 
isotropic,  v  is  Poisson’s  ratio  and  A  represents  the  magnitude  of  a  unit 
dislocation.  The  total  stress  at  t,o  is  then 
g(t)  =  t)  +  <7 

where  a  is  the  applied  stress  and  the  mean  value  of  g(t)  over  the  length  2r 
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2p  |  g{t)dt  =  g(a)  ^  g(a) 

Again,  in  the  presence  of  a  uniform  normal  tension  =  r  the  stresses 
developed  fay  a  pair  of  cylindrical  holes  radius  r  centered  at  x  =  ta  is  at 
points  x,o: 

r(x)  =  r  /  r2  +  3r4  +  r2  +  3r4  \  (12) 

2  |  (x-a) 4  (x-a)*  (x+a)2  (x+a)  *  J 

We  then  suppose  that  the  factor,  r  =  g(a),  is  known  and  in  the  first  instance 
determine  the  distribution  function  which  arises  from  the  action  of  the 
stress 

a  +  r(x) 

over  the  length  of  the  crack. 

Following  Head  and  Louat  [1]  we  have  the  distribution 

f(x)  =  -  1  1  l(e!-t!)V+r(t)  J  _dt  . 

^  J.c  t'x 

In  the  evaluation  of  this  integral,  it  is  helpful  to  note  that  we  can  write 
r(x)  =  r  f  x2  d  +  x4  d3  )  2a 

2  [  Ja  4’  3P  J  x*  -a2 

Substituting  for  r(x)  from  (12)  in  (13)  and  integrating  by  considering  the 
contour  integral 

(z2-t2)^  [c+r(z)]  dz_ 
z-x 


19 


we  find 


fix)  =  x 

7rA[c2-X2]^  I 


a  -  g(a)x2  fl  +  x2  d2  ) 
?  3a2 


f2  a  (2c2-a2-x2 

(a2-c2)^(a2-x2)2  , 


Then  according  to  (1)  the  stress  at  x  is  given  by 


=  1 


t r 


rc 


y-c 


x _  fa-  g(a)  (x2)  f  1  +  x2d2  ]  f  -  2a  (2c2-a2-x2)  )  dx  . 

-v2i F 1  L  l  (a2-c2)  ^  (a2-x2) 2  *‘a 


[c2*x2] 


Here  the  evaluation  of  the  contour  integral  is  complicated  by  the 
presence  of  a  triple  poles  at  z  =  a  and  a  double  pole  at  z  =  -a.  Neverthe¬ 
less,  the  algebra  is  straightforward  if  somewhat  lengthy  and  we  find  a 
total  stress 


a  =  a  +  g(a)  r2c2  f  1  +  r  d2  1  2a2-c2 

T"  l  wj  ^P)2 


But  by  definition  this  is  equal  to  g(a).  As  we  saw  above,  if  a-c/c  is 
small,  we  may  put  g(a).  On  this  basis  we  find 

g(a)  4a2(a2-c2)2  -  r2c2  (2a2-c2)  -  a  a 

{~rc7]i  • 


so  that 


g(a)  = 


era  4a2(a2-c2)2 

[  (4a2(a2-c ' 
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We  now  write  a  =  c  +  nr  where  nr/c<<i  and  obtain: 


g(a)  ^  era 

[a2-c2]’ 


l1  _re?7j 


and  have  from  (4)  and  with  r  =  g(a)r2 


f(x)  ~  XU 

rrA[c2-x2] 


7 


fl  -  r 

l  Wr- 


2a  2c2-a2-x2  .  a 
c7)-  (a2-x2)2  T-s 


}• 


where  s  =  _ 1_ 

i 


16n2 


The  evaluation  of  G  then  follows  from  Bilby  and  Eshelby’s  result  [5] 

G  =  ttAX  lim  (x-c)  f(x)2 
2  x-k: 


Substituting  we  obtain 


G  -  tto2c(1-v) 
-  TTCT2c(i-l') 


1  + 


1  1  V 

PI  J 


1  + 


ZI2  J 


so  that  G  is  increased  by  lOfc  when  n  =  2  and  the  cylinder  is  centered  at 
a  diameter  distance  from  the  tip  of  the  crack. 
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Surface  Displacement  and  Dynamic  Fracture 

Surface  displacement  near  the  tip  of  a  crack  has  become  of  importance 
with  the  use  of  the  so-called  caustic  method  of  examining  crack  behavior. 

In  this  method  the  lateral  surface  of  a  fracture  surface  is  polished  to  a 
mirror  finish.  The  area  adjoining  a  crack  tip  which  is  normal  to  this  sur¬ 
face  is  observed  using  a  system  of  optics  which  responds  to  changes  in 
inclination  of  the  surface.  Such  inclinations  are  induced  in  the  surface  as 
a  consequence  of  elastic  and  plastic  deformation  and  result  in  the  appear¬ 
ance  of  a  dark  area  bordering  the  crack  tip.  This  dark  area  or  rather  its 
boundary  is  the  caustic.  It  has  been  found  experimentally  that  this  curve 
approximates  to  a  circle.  This  description  is  found  to  be  increasingly  apt 
as  the  optics  are  adjusted  so  that  the  caustic  curve  is  generated  at  dis¬ 
tances  from  the  crack  tip  where  deformation  is  essentially  elastic. 

Again,  it  has  been  found  that  in  these  circumstances  that  there  is  a  one 
to  one  correspondence  between  the  crack  extension  force  applied  (as  deter¬ 
mined  from  Kjq)  and  the  radius  of  the  caustic.  This  allows  of  a  simple 
method  of  measuring  this  quantity. 

This  result  has  been  derived  theoretically  for  the  case  of  static  cracks 
and  also  for  a  dynamic  crack,  providing  its  speed  of  advance  is  signifi¬ 
cantly  less  than  that  of  sound. 

The  method  of  analysis  rests  on  the  representation  of  a  crack  in  mode  I 
as  an  array  of  edge  dislocations  in  the  plane  of  the  crack.  The  distribu¬ 
tion  of  this  array  is  heavily  concentrated  toward  the  tip.  Thus,  to  a  first 
approximation  the  crack  may  be  represented  as  ^  single  giant  dislocation 
located  at  the  crack  tip. 
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Now,  adopting  a  coordinate  system  in  which  the  axis  of  the  dislocation 
is  the  z-direction  while  the  plane  of  the  crack  is  represented  by  y  =  o,  we 
have  for  an  infinite  system  a  stress 


(7 


ZZ 


ubji+v)  y 
27t(1-v)  x^+y7 


Here  fj  is  the  shear  modulus,  v  is  Poisson’s  ratio  and  b  is  the  Burgers 
vector  of  the  dislocation.  When  such  a  dislocation  terminates  at  a  free 
surface  it  is  displaced  so  as  to  generate  stresses  which  just  cancel  those 
due  to  the  dislocation.  These  displacements  are  accompanied  by  inclina¬ 
tions  which  are  constant  along  lines  where  the  stress  is  constant.  It  is  a 
simple  matter  to  show  that  this  is  the  case  along  a  circle  whose  center 
lies  in  the  plane  of  the  crack  and  whose  circumference  passes  through  its 
tip. 

Thus,  theory  accords  with  experiment  in  confirming  the  usefulness  of 
the  caustic  approach. 
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PART  II  GRAIN  BOUNDARIES  AND  OSTWALD  RIPENING 

It  has  long  been  recognized  that,  other  things  being  equal,  materials  in 
which  grains  are  finer  are  tougher  than  those  where  the  grains  are 
coarser.  Grain  size  and  grain  growth  are  accordingly  germane  to  the 
study  of  fracture.  Furthermore,  the  rate  of  growth  of  grains  of  a  particu¬ 
lar  average  size  is  a  strong  function  of  the  size  of  particles  of  another 
phase  which  lie  in  the  grain  boundaries.  Broadly,  for  a  constant  volume 
fraction  the  restraint  offered  by  such  particles  decreases  as  the  size  of  the 
particles  increase  through  Ostwald  ripening.  Again,  there  is  evidence 
available  that  the  processes  of  normal  grain  growth  and  of  Ostwald  ripen¬ 
ing  (at  non-vanishing  concentration  of  second  phase)  are  both  stochastic  in 
origin. 

Accordingly  these  may  be  regarded  as  dual  topics.  Much  of  what  can  be 
said  for  one  is  applicable  to  the  other.  We  commence  with  grain  growth. 

Grain  Boundaries 

We  remark  that  in  a  polycrystal  not  ell  grains  are  of  the  same  size 
and,  in  fact,  grain  size  is  found  to  be  distributed  in  a  way  which,  remark¬ 
ably,  is  independent  of  mean  grain  size.  We  now  note  that  such  a  state  of 
affairs  could  be  expected  if  during  grain  growth  (or  Ostwald  ripening) 
grains  make  excursions  in  size  which  are  proportional  to  the  mean  and 
which  have  equal  probability  of  resulting  in  an  Increase  as  a  decrease  in 
size. 

On  this  basis  if  the  fraction  of  grains  of  size  lying  in  the  range  1  to 
1+<51  be  f  (1)  <51,  we  find  a  governing  equation 


af  _  a  a2f 
5t  5F 

where  A  is  constant. 
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(14a) 


For  this  an  appropriate  solution  is 

3(1)  =  1  1  e"1  /4At 

AF77* 

for  the  case  of  Ostwald  ripening  similar  considerations  hold.  These  will 
be  detailed  later. 

We  first  consider  how  this  approach  accords  with  the  whole  body  of 
evidence  relating  to  grain  growth. 

In  the  process  of  normal  grain  growth  some  grains  of  a  polycrystalline 
aggregate  vanish,  without  creating  voids.  The  volume  thus  made  available 
is  distributed  among  the  remaining  grains,  resulting  in  an  increase  in 
their  mean  size;  in  other  words,  in  grain  growth. 

Salient  features  of  this  phenomenon,  and  ones  with  which  any  accept¬ 
able  theory  must  be  in  accord,  are  that  the  processes  are  essentially 
steady  state  in  that  the  distribution  of  grain  sizes  is  invariant  with  the 
amount  of  growth  which  has  occurred,  that  the  rate  of  grain  growth  varies 
as  tn,  where  -l<r<-.5,  that  individual  grains  are  oriented  at  random,  and 
that  the  grain  structure  is  equi-axed. 

Measurements  relating  to  normal  grain  growth  have,  for  the  main  part, 
been  restricted  to  those  which  can  be  made  using  planar  sections.  In  par¬ 
ticular,  these  allow  the  determination  of  the  distribution  of  grain  sizes 
and  of  the  number  of  faces  to  a  grain,  the  rate  of  grain  growth  and  its 
dependence  on  time.  Similar  results  have  also  been  obtained  from 
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computer  modeling.  It  is  required  of  theory  that  it  give  correct  qualitive 
predictions  of  the  results  and  also  those  obtained  when  the  measurements 
relate  to  three  dimensions. 

A  first  task  in  delineating  the  requirements  on  theory  is  to  decide  just 
what  is  the  mathematical  form  of  the  distribution  of  grain  size,  1.  In  the 
past,  this  distribution  has  been  fitted  to  four  different  functions:  lognor¬ 
mal^,  the  Rayleigh  function  [8,9],  a  complex  exponential  forrrJ^  and 
that  of  Lifshitz  and  Slyozov^  These  functions  are  respectively: 


g(l) 

_  A  e-ln3(l/lo)/2<T*  . 

(15) 

f(l) 

=  B  1  .  ; 

(16) 

and 

h(l) 

=  ce-D([l]J-d>*  . 

k(l) 

=  /9u  .  (2e/  e  '2^/2'u  ,  u  =  \/\0 

TK  *19  -  A 

(17) 

u)  *  + 


fi  is  equal  to  the  dimensionality  involved,  i.e.  2  or  3. 

Here  A,  B,  C,  D,  10,  and  1,  are  disposable  constants. 

No  one  of  these  equations  gives  a  good  fit  to  all  the  data.  The  best 

individual  fit  is  probably  that  of  (17)  applied  to  the  data  given  by  Aboav 
f  131 

and  Langdon  1  Again,  it  would  seem  that  this  equation  is  capable, 
through  adjustments  to  the  parameter  D,  of  giving  a  reasonably  acceptable 
fit  to  all  the  data.  However,  since  any  theoretical  justification  seems 
improbable,  it  is  difficult  to  accept  that  this  agreement  is  other  than 
fortuitous.  The  difficulty  with  this  form  is  essential  and  lies  in  the  fact 
that  the  distribution  does  not  tend  to  zero  at  zero  grain  size 

In  recent  years  there  has  been  a  growing  conviction  among  many 
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workers  that  the  distribution  of  grain  sizes  approximates  to  the  lognormal 
form.  This  conclusion  usually  rests  on  the  following  observations.  Data 
is  plotted  using  special  graph  paper  in  which  the  abscissa  is  the  cumula¬ 
tive  percentage  of  grains  measured.  This  scale  is  non-linear  and  given  by 


F(x)  = 


£-(ln  d/l0)  V2c72 


dl/1 


J  x 


“-(In  l/i0)2/2cr2 
.  o 


dl 

T 


The  ordinate  is  scaled  as  lnx.  Data  which  obeys  a  lognormal  distribution 
then  plots  as  a  straight  line  with  slope,  o. 

Thus,  an  examination  of  such  plots  shows  a  small  but  definite  deviation 
from  the  straight  line  configuration.  Data  are,  as  is  frequently  recog¬ 
nized,  only  approximately  lognormal.  Consistently,  data  plotted  directly 
against  the  logarithm  of  the  grain  size  shows  significant  assymetry  about 
the  peak.  Again,  a  direct  comparison  of  the  relations  (2)  and  (3)  shows 
that  they  can  be  In  essential  accord  if  the  constants  B,  C,  D  and  have 
suitable  values.  It  was  shown  that  the  two  are  in  reasonable  agreement 
when  a  .5.  This  suggests  that  a  plot  of  data  which  actually  follow  the 
Rayleigh  distribution  would  approximate  to  a  straight  line  when  plotted  on 
the  special  graph  paper  mentioned  above.  Clearly,  the  data  of  Conrad, 
Swlntowski  and  Mannan  1  J  do  conform  approximate'^  to  straight  lines 
which  do  not  deviate  far  from  the  curves.  Nevertheless,  the  latter 
provide  a  better  fit  to  the  data  than  do  the  lines.  The  standard  deviation 
in  the  lognormal  representation  of  these  results  approximates  to  .6  and 
it  is  apparent  that  this  is  the  case  for  most  published  material.  Thus, 
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Schukler^^  lists  results  from  eight  materials.  One  of  these  cr-iron,  is 
described  as  not  not  being  lognormal,  the  remaining  seven  have  standard 
deviations  which  average  .59  and  which  lie  in  the  range  0.468  to  0.674. 
We  conclude  that  all  experimental  results  which  give  standard  deviations 
of  about  this  amount  on  the  lognormal  approximation  will  be  in  reasonable 
agreement  with  the  Rayleigh  distribution.  This  is  not  to  say  that  all 
grain  size  distributions  describable  as  lognormal  can  also  be  considered 
to  be  Rayleigh.  For  example,  Rhines  and  Patterson  have  found 
distributions  in  which  a  -  2.  Data  approximating  to  a  lognormal  distri¬ 
bution  for  standard  deviations  in  this  range  is  not  simply  describable  in 
terms  of  Rayleigh.  However,  it  should  be  borne  in  mind  that  these  data 
cannot  be  expected  to  be  representative  of  normal  grain  growth  since  they 
were  obtained  using  pre-strains  in  the  range  2  -  30*;  values  so  low  as  to 
give  non-uniform  nucleation  during  recrystallization. 

Passing  now  to  a  consideration  of  grain  shape,  that  is  the  distribution 
of  a  number  of  faces,  or  in  two  dimensions,  of  edges  per  grain,  it  is  clear 
that  here  the  form  is  lognormal  or  a  very  good  approximation  to  it.  All 
authors,  including  Hu,  who  found  that  his  grain  size  distribution  data  was 
not  lognormal,  appear  to  agree  on  this  point.  This  discrepancy  is  not  at 
variance  with  the  several  relations  which  have  been  found  between  the 
diameter  of  a  grain  and  the  number  of  its  sides: 
n  =  3  +  kl 

where  n  is  the  average  number  of  sides  of  per  grain  of  diameter  1^; 
and 


N  =  2  +  KL 
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where  L  is  the  average  diameter  of  a  grain  with  an  integral  number  of 
sides  If  all  grains  having  N  sides  were  the  same  size,  it  would  be  ex¬ 
pected  that  grain  size  and  shape  would  have  the  same  distribution.  As  it 
is,  they  should  be  similar  but  not  identical. 

Allowing  that  there  is  this  difference,  no  explanation  has  yet  been  ad¬ 
vanced  as  to  why  the  distribution  of  shape  should  be  lognormal.  Toward 
such  an  explanation  we  recognize  that  observations  of  grain  shape  derive 
from  planar  intersections  of  three-dimensional  ensembles.  In  contra¬ 
distinction  with  grain  size  there  is  no  unique  transform  relating  charac¬ 
teristics  in  two  to  those  in  three  dimensions.  Thus,  the  shape  observed 
as  a  consequence  of  a  planar  intersection  with  one  particular  grain 
depends  not  only  on  the  number  of  its  faces  but  also  on  their  sizes  and 
orientation  with  respect  to  the  intersecting  plane;  that  is  to  say,  a  multi¬ 
plicity  of  transforms  are  operative.  Then  allowing,  as  has  been  shown^, 
that  if  the  distribution  of  grain  size  is  Rayleigh  in  two  dimensions  it  is 
Rayleigh  in  three  also,  we  see  that  the  distribution  of  shape  in  two- 
dimensions  cannot  be  Rayleigh  if,  as  seems  reasonable,  we  suppose  it  is 
in  three.  We  now  suppose  that  the  observed  shape  distribution  is  a  linear 
superposition  of  different  and  independent  distributions  each  derived  from 
the  operation  of  a  different  transform  operating  on  the  distribution  of 
shape  in  three  dimensions.  Each  of  these  distributions  must  accord  with 
the  Euler  rules  so  that  the  average  number  of  sides  to  a  grain  must  be  6. 
Again,  supposing  that  no  grain  have  less  than  three  sides,  we  require  that 
each  distribution  vanish  when  n,  the  number  of  sides,  is  2  and  when  it  is 
infinity. 
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Now,  according  to  the  central  limit  theorem  of  probability  theory  these 
characteristics  are  consistent  with  the  appearance  of  a  lognormal  distribu¬ 
tion  and  it  is  to  this  that  we  address  its  appearance. 

The  Rayleigh  distribution  was  derived  on  the  basis  that  grain  growth  is 
a  stochastic  process.  Besides  the  accord  between  theory  and  experiment 
indicated  above,  additional  support  for  this  approach  has  recently  become 
available  as  the  result  of  computer  modeling  of  grain  growth.  This  sup¬ 
port  comes  in  four  ways: 

The  calculated  distribution  is  in  good  agreement  with  the  Rayleigh 
nearly  everywhere;  the  distribution  is  found  to  be  the  same  in  three  as  it 
is  in  two  dimensions.  This  is  critical*  since  it  is  a  unique  property  of  the 
Rayleigh  function  that  it  is  invariant  to  the  Abelian  transform  used  to 
generate  the  distribution  function  in  three  dimensions  from  that  in  two,  or 
vice  verse; 

Consideration  of  the  trajectories  of  individual  grains  in  grain  size 
space  shows  that  size  and  growth  rate  are  not  directly  related.  Further^ 
more,  as  postulated,  grains  do  show  changes  in  the  direction  of  their 
growth  even  though  they  are  much  larger  than  the  mean  grain  size. 

Finally,  this  work  also  indicates  a  linear  relation  between  the  number 
N  of  sides  of  a  grain  and  the  mean  diameter  L  of  grains  of  this  class. 

Thus, 

N  =  2  +  KL 

The  constant,  K  was  evaluated  and  found  to  be  0.23.  Assuming  that  such  a 


*  Other  functions  will,  however,  satisfy  this  transform  to  various  degrees 
of  approximation. 
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relation  exists,  the  value  of  K  can  be  evaluated  from  the  Rayleigh  distri¬ 
bution  on  the  basis  that  the  average  corresponds  to  grains  with  six 
fl71 

sides1  .  The  value  predicted  in  this  way  is  0.25. 

Thus,  there  is  now  a  significant  amount  of  evidence  to  support  the  view 
that  random  walk  is  of  dominating  importance  in  grain  growth.  According¬ 
ly,  we  shall  examine  its  basic  features  so  as  to  determine  whether  such 
behavior  should  be  anticipated.  We  shall  then  be  concerned  to  justify  the 
couchment  of  resulting  analysis  in  terms  of  the  linear  dimensions  of  the 
grains  rather  than,  say,  their  volume. 

In  this  connection,  the  first  matter  to  be  settled  is  whether  grain 
growth  is  to  be  regarded  as  involving  the  grains  in  and  of  themselves  or 
simply  the  interfaces  between  them.  We  are  inclined  to  the  view  that 
since  only  the  boundaries  move  during  grain  growth  the  presence  or 
absence  of  the  material  within  the  grain,  which  makes  up  the  grains,  is 
largely,  if  not  entirely,  irrelevant.  Thus,  we  expect  the  constraint  of 
constant  volume  (area  in  two  dimensions)  while  important  and  necessarily 
included,  need  not  be  central  to  the  theory.  Rather,  we  should  include 
those  constraints  which  arise  from  the  fact  that  we  are  concerned  with  an 
assembly  of  faces  in  quasi -equilibrium  under  the  forces  due  to  their 
intrinsic  energy;  that  is  to  say,  we  take  due  note  of  the  so-called  Euler 

nW17l. 

In  three  dimensions,  these  require  that  3  faces  meet  in  a  line  and  four 
such  lines  meet  at  a  point  forming  a  four-fold  node.  In  two  dimensions, 
three  edges  must  meet  at  a  point.  On  this  basis  and  with  the  additional 
constraint  that  edges  at  a  node  are  spaced  120°  apart,  Von  Neuman^ 
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examined  the  behavior  of  an  n-cornered  cell  of  a  two-dimensional  bubble 

froth.  From  his  determination  of  the  resultant  forces,  it  can  be 
r  1 2 

inferred1  ]  that  the  cell  will  tend  to  grow  at  a  rate 

dl  _  M  n-6  (18) 

dt  n 

where  M  is  a  constant,  n  the  number  of  comers  (or  sides)  and  1  a  measure 
of  the  linear  dimensions  of  the  cell.  A  comparable  expression  applicable 
to  three  dimensions  has  yet  to  be  determined,  but  it  can  be  expected  that  a 
somewhat  similar  equation  should  hold  there,  perhaps 

dl  _  K  m-24  (19) 

3t  m 

where  m  is  the  number  of  grain  edges.  The  numbers  6  and  24  of  eqs.  (18) 

and  (19)  are  analogous.  In  two  dimensions,  an  array  of  grains  obeying 

modified  Euler  rules  is  stable  when  each  grain  has  six  comers^^.  In 

n  9] 

three  dimensions  a  stable  space-filling  ensemble  of  polyhedra  is  found'  J 
when  each  has  14  faces,  24  edges  and  36  apices.  Thus,  much  if  not  all 
the  physics  of  grain  boundary  structure  essential  to  our  purposes  is  encap¬ 
sulated  in  (18)  and  (19)  and  it  is  important  to  stress  that  these  equations 
relate  numbers  specifying  a  topological  quantity  to  a  linear  dimension. 
Allowing  that  grain  growth  involves  an  increase  in  size  of  the  average 
grain  of  an  assembly,  it  is  apparent  that  such  growth  can  occur  only  as  a 
result  of  the  disappearance  of  particular  grains.  We  can  as  direct  conse¬ 
quence  of  the  modified  Euler  conditions  identify  these  disappearing  grains 
as  tetrahedra  (triangles  in  two  dimensions).  Thus,  a  tetrahedron  is 
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unique  among  polyhedra  in  having  four  apices.  Migration  of  these  apices 
to  a  central  point  leaves,  as  required  for  stability,  a  four-fold  node.  This 
fact,  together  with  (19)  would  seem  to  suggest  a  model  in  which  a  grain 
with  fewer  than  24  edges  (six  comers)  would  progressively  lose  them 
until  it  finally  disappears,  while  conversely,  grains  with  more  than  this 
number  grow  indefinitely.  These  features  together  with  the  seemingly 
reasonable  assumption  that  grain  diameter  is  proportional  to  n  forms  the 
basis  for  the  treatment  given  by  HillerJ^  and  expressed  by  either  of  the 
eouations: 

t 


where  a  and  b  are  constants.  1  is  the  diameter  of  a  six  sided  grain. 

An  immediate  difficulty  with  this  approach  is  that  it  takes  no  account 
of  grain  contiguity;  an  omission  which  we  now  attempt  to  rectify.  In  the 
interests  of  simplicity  we  restrict  the  discussion  to  considerations  in  two 
dimensions.  We  observe  that  while  von  Neumann’s  analysis  indicates 
whether  a  particular  grain  will  tend  to  grow  or  to  shrink,  It  does  not  bear 
on  the  question  as  to  whether  or  not  a  particular  side  will  extend  or  con¬ 
tract,  possibly  vanishing.  This  possibility  arises  because  each  side  is 
(1)  shared  by  two  grains  and  (2)  terminates  in  two  grains.  Toward  a 
determination  of  the  consequences  of  these  fee '-res.  we  first  consider  the 
well  known  change  in  configuration  illustrated  in  Fig.  2. 
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FIG.  2 

Illustrating  a  change  in  the  grain  configuration  in  two  dimensions. 

Here  the  merger  of  the  grain  comers  x  and  y  and  their  replacement  by  p 
and  t  is  accompanied  by  the  loss  of  one  of  one  side,  by  each  of  the  grains  A 
and  B  and  an  equal  gain  by  both  C  and  D.  The  converse  is  true  if  p  and  t 
merge.  As  demanded  by  Euler’s  rule  there  is  not  change  in  the  number  of 
sides  in  either  case.  Thus,  a  grain,  for  example  A,  may  either  gain  or 
lose  a  side  by  one  or  other  of  the  procedures.  A  question,  crucial  in  the 
assumption  expressed  in  (6)  now  arises:  to  what  extent  is  the  choice 
between  these  alternatives  determined  by  the  characteristics  of  grain  A? 
Toward  a  resolution  of  this  question,  we  first  examine  the  average  be¬ 
havior  of  an  n-sided  grain  and  refer  to  the  grain  comer  illustrated  in 
Fig.  2.  Here  the  lines  represent  the  chords  joining  the  ends  of  grain  edges 
which  are  taken  to  be  bent  into  arcs  of  circles 
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FIG.  3 

Grain  Corner  Geometry 


We  consider  the  cases  where  n  >  6  so  that  a ,  the  mean  value  a,  is  greater 
than  77-/6.  We  see  that  as  shown,  the  point  O  is  not  at  a  position  of 
equilibrium  and  consequently  that  it  will  move  in  response  to  the  line 
tensions  T  directed  along  the  directions  OP,  etc.  as  shown.  If  the  result¬ 
ant  motion  in  the  direction  OS  lies  within  the  right  angle  QOR,  the  grain 
edge  will  tend  to  shorten;  if  it  lies  within  the  angle  ROT,  the  edge  will 
lengthen.  Resolving  forces,  the  condition  for  shortening  is  easily  shown 
to  be 

T  sin  y _  >  cos  tor. 

T  cos  y  -  2T  cos  or 

When  rearranged  this  becomes 


cos(y  +  or)  <2  cos2a. 


(21) 
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This  condition  is  satisfied  except  when  y  may  be  regarded  as  small.  Thus, 
when  n  =  7,  a  =  64.29  we  find  that  on  average  the  side  OQ  tends  to 
shorten  unless  y<yc,  yc  =  3.59. 

The  most  restrictive  case  occurs  when  y_  =  y  ,  the  grain  in  POQ  is 

c  max 

12-sided  and  yc  =  7.5°.  Similar  considerations  can  be  applied  to  the  case 
where  grains  have  less  than  six  sides.  There,  we  find,  when  a  takes  its 
mean  value,  a,  the  sides  OQ  move  so  as  to  shorten,  for  all  y. 

It  is  helpful  in  appreciating  the  implications  of  these  results  to  con¬ 
sider  the  case  of  a  six-sided  grain  in  the  special  circumstances  where  each 
angle  POQ  =  2a  has  the  mean  value  2n/3.  Here  (21)  is  satisfied  for  all 
positive  y  and  fails  when  y  is  negative.  But  at  each  grain  comer,  there 
are  two  grain  sides  to  be  considered.  In  these  two  cases  y  takes  equal  and 
opposite  values.  Accordingly,  one  side  shrinks,  the  other  lengthens.  Then 
since  each  edge  has  two  ends  we  see  that  four  possible  situations  can 
arise.  In  one,  both  ends  tend  to  move  together,  in  another  they  move 
apart,  in  the  remainder  the  motions  oppose.  If  then,  as  seems  to  be 
indicated  we  assign  equal  probability  to  each  case  we  see  that  the  chance 
of  a  side  lengthening  is  just  equal  to  that  of  its  shortening.  Thus,  as 
might  be  expected,  we  find  that  a  six-sided  grain  has  equal  probility  of 
tending  to  gain  as  to  lose  sides. 

In  the  case  of  grain  with  more  than  six  sides  there  is  bias  toward  the 
extension  of  all  sides.  The  amount  of  this  bias  is 

tc/6  = R 

8  is  the  range  of  6  (see  Fig.  2)  found  from  an  appropriate  weighting  for 


the  frequency  of  occurence  of  particular  values.  Here  it  is  sufficient  to 
recognize  that  0  can  lie  in  the  range  60°  —  150°  and  hence,  in  view  of 
(21),  that  R  should  be  significantly  less  than  1  in  all  cases  in  which  n>6. 
We  conclude,  then,  that  the  processes  by  which  grains  gain  and  lose  sides 
during  growth  are  in  the  main  random  rather  than  directed  when  n>6. 

When,  on  the  other  hand,  n<6  the  situation  is  reversed.  With  decreas¬ 
ing  n  there  is  then  a  steadily  increasing  correlation  between  a  change  in 
grain  size  and  the  number  of  its  sides.  This  is  so  because  (21)  is  satis¬ 
fied  for  all  y  when  a  takes  its  mean  value  a,  since  cKtt/ 3.  Thus,  (21) 
can  fail  and  the  grains  have  a  tendency  to  gain  a  side  only  at  those  corners 
at  which  stochastic  variation  provides  a  value  of  cr  exceeding  tt/3.  Such 
excursions  can  be  expected  at  reasonably  common  when  n  =  5  and  or  =  w/2 
but  rare  when  n  =  A  (a  =  7r/6).  However,  Imam,‘^'  using  a  p'noto- 
emission  electron  microscope  to  observe  grain  boundary  motion  in  titan¬ 
ium  has  reported  a  sequence  in  which  a  shrinking  four-sided  grain  gained 
a  side  before  losing  it  again. 

Thus,  we  find  that  when  n>6,  the  path  of  addition  (and  subtraction)  of 
grain  sides  should  be  essentially  that  of  a  random  walker.  The  dominance 
of  directed  "flow"  by  stochastic  factors  is  continuously  reduced  and  finally 
disappears  as  n  decreases  below  6. 

The  effect  of  this  variation  has  been  examined  by  Pande1  -  who  has 
modeled  this  system  by  considering  the  equation 

3f  _  A  32f  +  B  3(f/l)  . 

3t  cfF  "31 — 

Here  the  second  term  represents  the  effect  of  directed  flow.  This  equation 
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has  a  closed  form  solution.  For  reasonable  values  of  the  ratio  3/A  the 
distribution  deviates  only  slightly  from  that  calculated  when  B  =  0,  namely 
that  of  Ravieigh. 

We  see  then  that  the  Rayleigh  distribution  is  to  be  exoected  on  the  basis 
of  the  random- walker- view  of  grain  growth  and  that  is  observed  exoeri- 
mentaily.  We  conclude,  then,  from  this  agreement  and  other  associated 
accords  that  there  is  strong  evidence  to  supoort  that  the  mechanism  asso¬ 
ciated  with  this  view  is  indeed  operative  and  controlling  in  grain  growth. 

Ostwald  Ripening 

In  place  of  'Ida)  we  have,  again  a  basis  of  random  excursions: 

m  r  ■) 

24  =  B  Ax  "  Srd. 

where  x  =  dtrc2,  r  is  the  radius  of  a  particle,  assumed  spherical  and  B  is  a 
rate  factor,  fix'  is  now  the  distribution  of  surface  areas  of  the  particles. 
Surface  area  rather  than  particle  radius  is  adopted  (without  prejudice  to 
the  result)  as  being  the  more  physically  realistic  quantity. 

To  fix  3  we  recognize  that  surface  displacements  scaie  with  the  size  of 
the  surface  and  thus  with  r:  while  the  rate  of  change  of  particle  radius 
varies  as  1  /r. 

Thus . 


3  ~ 


where  r  and  x  are  the  mean  values  of  r  and  x,  and  (22)  becomes: 


(22) 
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3f  _  A  x*  32f 


St  " 


3  *} 


where  A  is  a  constant. 


Then  supposing  x  ^  x(t)  and  that 


rt 


r  = 


x^(t)  dt 


(23)  becomes 


3f 

Si 


A  3ff 
3x2 


with  a  solution 


f  =  x 


-  x 


Then,  rate  of  growth  of  x  is  as  r*  and 


dr  _  S*  orr1/4,  r3/4a  t,  r^a  t1/3 

3T 


(23) 


(24) 


(25) 


Thus  we  conclude  on  this  basis  that  Ostwald  ripening  proceeds  as  the  one- 
third  power  of  the  true.  This  accords  with  experiment. 

Again,  expressed  in  ordinary  time  units  and  particle  radius,  the 
distribution  becomes 

f  =  rf  e  '  ( t™  )  /4A 
t 


This  distribution  has  been  found  to  accord  well  with  available  experimental 
evidence. 
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